Lepton flavour violating processes at the International Linear Collider 



P.M. Ferreira 1 *, R.B. Gucdes 1 t an d R. Santos 1 ' 2 * 
1 Centra de Fisica Teorica e Cornputacional, 
Faculdade de Ciencias, Universidade de Lisboa, 
Avenida Professor Gama Pinto, 2, 1649-003 Lisboa, Portugal 
2 Department of Physics, Royal Holloway, University of London, 
Egham, Surrey TW20 OEX United Kingdom 

(Dated: November, 2006) 

We study the effects of dimension six effective operators on the flavour violating production and 
decay of leptons at the International Linear Collider. Analytic expressions for the cross sections, 
decay widths and asymmetries of all flavour changing processes will be presented, as well as an 
' analysis of the feasibility of their observation at the ILC. 
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I. INTRODUCTION 



In the next years the Large Hadron Collider (LHC) will start to function and provide the scientific 
community with a new tool with which to explore hitherto unknown regions of particle physics. We expect 
many exciting discoveries to arise from LHC experiments. However, the LHC is a hadronic machine, and 
as such precision measurements will be quite hard to undertake there. Also, the existence of immense 
backgrounds at the LHC may hinder discoveries of new physical phenomena already possible at the 
■ energies that this accelerator will achieve. Thus it has been proposed to build a new electron-positron 

collider, the International Linear Collider This would be a collider with energies on the TeV range, 
with extremely high luminosities. The potential for new physics with such a machine is immense. In 
this paper, we will focus on a specific sector: the possibility of processes which violate lepton flavour 
occurring. 

We now know that the solar and atmospheric neutrino problems 0] arise, not from shortcomings of 
solar models, but from particle physics. Namely, the recent findings by the SNO collaboration Q have 
shown beyond doubt that neutrinos oscillate between families as they propagate over long distances. 
Leptonic flavour violation (LFV) is therefore an established experimental fact. The simplest explanation 
for neutrino oscillations is that neutrinos have masses different from zero - extremely low masses, but 
non-zero nonetheless. Oscillations with zero neutrino masses are possible, but only in esoteric models Q. 
With non-zero neutrino masses, flavour violation in the charged leptonic sector becomes a reality (whereas 
with massless neutrinos, it is not allowed in the Standard Model (SM)). This is a sector of particle physics 
for which we already have many experimental results [j| , which set stringent limits on the extent of flavour 
violation that may occur. Nevertheless, as we will show in this paper, even with all known experimental 
constraints it is possible that signals of LFV may be observed at the ILC, taking advantage of the large 
luminosities planned for that machine. There has been much attention devoted to this subject. For 
instance, in refs. Q effective operators were used to describe LFV decays of the Z boson. LFV decays 
of the Z boson were also studied in many extensions of the SM [7|. The authors of refs. [alSl performed 
a detailed study of LFV at future linear colliders, originating from Supersymmetric models. Finally, a 
detailed study of the four fermion operators in the framework of LFV is performed in 10]. In that work 
the exact number of independent four fermion operators is determined. Gauge invariance is then used to 
constrain LFV processes which are poorly measured, or not measured at all. 

In this work we carry out a model-independent analysis of all possible LFV interactions which might 
arise in extensions of the SM. To do so, we utilise the effective operator formalism of Buchmuller and 
Wyler [Tlj ] . a standard tool in such studies. This formalism parameterizes whatever new physics may 
appear in theories that generalise the SM as effective operators of dimension greater than four. Our 
goal is to provide the reader with as many analytical expressions as possible for physical quantities 
which might be measured at the ILC. In this way, our experimentalist colleagues will have expressions 
they can include in their Monte Carlo simulations. This paper is structured as follows. In section [TT] 
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we present the effective operator formalism and list the operators which contribute to lepton-violating 
interactions, both interactions with gauge bosons and four-fermion contact terms. In section Hill we use 
the existing experimental bounds on decays such as /i — ► e 7 to exclude several of the operators which 
could a priori have a contribution to the processes we will be considering. We also analyse the role that 
the equations of motion of the fields play in further simplifying our calculations. Having chosen a set 
of effective operators, we proceed, in section IIV1 to calculate their impact on LFV decays of leptons, 
deducing analytical expressions for those quantities. Likewise, in section [Vl we will present analytical 
results for the cross sections and asymmetries of several LFV processes which might occur at the ILC. 
We analyse these results in section I VII performing a scan of a wide range of values for the anomalous 
couplings we introduced, and considering their possible observability at the ILC. 



II. FLAVOUR CHANGING EFFECTIVE OPERATORS 



The effective operator formalism of Buchmiiller and Wyler [ll| is based on the assumption that the 
Standard Model (SM) of particle physics is the low energy limit of a more general theory. Such theory 
would be valid at very high energies but, at a lower energy scale A, we would only perceive its effects 
through a set of effective operators of dimensions higher than four. Those operators would obey the 
gauge symmetries of the SM, and be suppressed by powers of A. This allows us to write the effective 
lagrangian as a series, such that 

C = C*» + ±cV + jsCM+o(J^ , (1) 

where L SM is the SM lagrangian and and contain all the dimension five and six operators which, 
like C SM , are invariant under the gauge symmetries of the SM. The terms break baryon and lepton 
numbers. Hence, we should start by considering the dimension five LFV terms. However, these terms 
would also be responsible for the generation of neutrino masses. With the present limits on neutrino 
masses, m„ < 2 eV @, the scale of new physics would have to be of the order 10 13 GeV 11 1, which is 



clearly out of the reach of the next colliders. This leaves us with the operators, some of which, after 
spontaneous symmetry breaking, generate dimension five terms. The list of dimension six operators is 
quite vast [ll|. In this work we are interested in those operators that give rise to LFV. Throughout this 
paper we will use lh to represent a heavy lepton and k denotes a light one (whose mass we consider zero) . 
In processes where a tau lepton is present, both the muon and the electron will be taken to be massless. 
If a given process only involves muons and electrons, then the electron mass will be set to zero, but the 
muon mass will be kept. Whenever the lepton's mass has no bearing on the result we will use I for all 
massless leptons, and drop the generation index. 

The effective operators that will be important for our studies fall in three categories: (a) those that 
generate flavour- violating vertices of the form Zlhh and jlhh (and also, for some operators, vertices 
like "/"flhh)', these operators always involve gauge fields, either explicitly or in the form of covariant 
derivatives, (b) Four-fermion operators, involving only leptonic spinors. (c) And a type of operator that 
involves only scalar and fermionic fields that will roughly correspond to a wave function renormalization 
of the fermion fields. 



A. Effective operators generating Zlhh and jlhll vertices 



There are five tree-level dimension 6 effective operators that can generate a new Z lh k interaction. This 
means that these interactions are compatible with SM symmetries at tree level. Following the notations 
of [ll| we write the first two operators as 

Od c = (4^ M 4) , De = || [dH\ 4) . (2) 

The coefficients 77^ are complex dimensionless couplings and the (i, j) are flavour indices. For flavour 
violation to occur, these indices must differ. t L is a left-handed SU(2) doublet, e 3 R is a right-handed 
U(l)y singlet, 4> is the Higgs scalar SU{2) doublet. Notice that the terms contributing to the interaction 
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Z lh k in which we are interested appear in the lagrangian when the Higgs doublet acquires a vacuum 
expectation value (vev) v. There is no jlhh interaction stemming from these terms, although one may 
obtain contributions to vertices involving also a Higgs field, such as 7 <plh h and Z <plhh- 

The remaining three operators that contribute to the vertices Z lh h but not to 7, lh h are given by 

e L{1) / \ e L{3) /_ x 

°fz= * -JT (^M(47 M 4) . Og>= i -^5- (^2W) (4tV4) . (3) 

Again, ^ and 0^- are complex dimensionless couplings, and the contributions to Z lh h arise when 
both scalar fields acquire a vev v. Because the covariant derivatives act on those same fields and the SM 
Higgs has no coupling to the photon, there are no contributions to 7 lh h from these operators. There are 
however five dimension six operators that contribute to both the Z lh h and 7 lh h vertices and are only 
present at the one-loop level. They are given by 

a BR 



a BL /_ N a WL / . \ 

o f .B = 1^ (i l L Y L d v e L ) , o m = i (t L rj r d v e L ) wl 

aW 

(ll 4) B, v , O eW<t> = (f L TI a** 4) <f> . (4) 



O 



eB<p 



ctij and Pij are complex dimensionless couplings, B^ IV and Wl v are the usual U(l)y and SU(2)l field 
tensors, respectively. These tensors "contain" both the photon and Z boson fields, through the well- 
known Weinberg rotation. Thus they contribute to both Z lh k and 7 lh h when we consider the partial 
derivative of in the equations (j4|) or when we replace the Higgs field cf> by its vev v in them. We will 
return to this point in section fin Al 



B. Four-Fermion effective operators producing an eelhh contact interaction 



Because we are specifically interested in studying the phenomenology of the 1LC, we will only consider 
four-fermion operators where two of the spinors involved correspond to the colliding electrons/positrons 
of that collider. Another spinor will correspond to a heavy lepton, lh- There are four relevant types of 
four-fermion operators that contribute to e + e~ — > lh h, 

Oft = "-f (h^ L ) («) , Cg* = K \ {hl^lL) (h*fT*t L ) , 

O ee = ~y (eRj^e R ) (e R ^e R ) , O le = n le (l L e R ) (e R £ L ) . (5) 

Again, all of the couplings in these operators are, in general, complex. As we have done with the previous 
operators, we should now consider all possible "placements" of the lh spinor, and consider different 
couplings for each of them. But that would lead to an unmanageable number of fermionic operators, 
all with the same Lorentz structure but differing simply in the location of the heavy lepton spinor. 
Thus we will simplify our approach and define only one coupling constant for each type of operator. 
An exception is the operator Oe e = Kg e (IlCr) (Sr£l), which corresponds to an interaction between a 
right-handed current and a left-handed one. Depending on where we place the lh spinor, then, we might 
have two different effective operators. For example, if we consider the operators that would contribute 
to e + e~ — > r~ e + , the two possibilities we would have, putting the chiral structure of the operators in 
evidence, are 

Te = K Te (t L ^ R e R ) (e R ^ L £ e R ) , O eT = K eT (1% 7fl 6 R ) (f R f L £%) , (6) 

where £ e and £ T are the leptonic doublets from the first and third generations, respectively, and "/l,r = 
(1 T 7s)/2 are the usual chiral projectors. As we see, we find two different Lorentz structures depending 
on where we "insert" the r spinor. Therefore we define two different couplings, each corresponding to 
the two possible flavour- violating interactions. 
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It will be simpler, however, toparameterize the four-fermion effective Lagrangian built with the operators 
above in the manner of ref. [13. For the e + e~lhh interaction, we have 



C, 



I,J=L,R 



Vfj (e 7M 7/e) (JhTfuk) + Sfj (e 7/ e) (kjjh) 



(7) 



The vector-like (Vj j) and scalar-like (SVj) couplings may be expressed in terms of the coefficients of the 
four four-fermion operators written in eq. ^ [13j in the following manner: 



LL 



- (k {1) 

2 r« 



.(3) 



Vrr 



V L r = 0, V rl = 0, 



Srr = , Sll — , Slr. — Kg e , S_ 



RL 



(8) 



C. Effective operators generating an lh l\ mixing 



There is a special kind of interaction that corresponds to a wave-function renormalization, which has its 
origin in the operator 

O eH = , (9) 

where Sij are complex dimensionless couplings. After spontaneous symmetry breaking the neutral compo- 
nent of the field <f> acquires a vev (</>o — > </>o + v, with v — 246/ \/2 GeV) and a dimension three operator 
is generated which is a flavour-violating self-energy like term. In other words, it mixes, at the level of 
the propagator, the leptons of different families. We consider these operators here for completeness, even 
though we will show that they have no impact in the phenomenology whatsoever. 



III. THE COMPLETE LAGRANGIAN 

The complete effective lagrangian can now be written as a function of the operators defined in the previous 
section 

c = c eellh + o De + o Bc + + o$ + of] + 

+ O eB + O lB + O ew + O eB< t> + O eW 4> + h.c. . (10) 

This lagrangian describes new vertices of the form jlhh, Zl^h, eelhh, hih (and many others) and 
all of their charge conjugate vertices. We will also consider an analogous lagrangian with flavour indices 
exchanged - in other words, we will consider couplings of the form rjhi and rjih, for instance - except for the 
four fermion lagrangian, as was explained in the previous section. Rather than write the Feynman rules 
for these anomalous vertices and start the calculation of all LFV decay widths and cross sections, we shall 
use all experimental and theoretical constraints to reduce as much as possible the number of independent 
couplings. After imposing these constraints we will write the Feynman rules for the remaining lagrangian 
and proceed with the calculation. 



A. The constraints from lh — * h 7 



Some of the operators presented in the previous section can be immediately discarded due to the very 
stringent experimental bounds which exist for the decays t —> [ij, t — > and fi — > ej. The argument 
is as follows: all the operators in eqs. ([4J contribute to both 7 lh li and Z lh h interactions, due to the 
presence of the gauge fields and in the field tensors -B M „ and W^ v that compose them. Then we 
can write, for instance, an operator C^ 7 , given by 

O tj = (t L rd v £{) (11) 
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where is the usual electromagnetic tensor. This operator was constructed from both (Dig and Oew, 
and the new effective coupling a]j L is related to a B L and L through the Weinberg angle 8w by 



7 L 

I ■ ■ = 



>3 



(12) 



Following the same exact procedure we can also obtain an operator Ogz-, with coupling constant given 

by 



a 



Z L 



sm vw ol 



BL 



COS O \y a 



W L 

ij 



(13) 



New operators with photon and Z interactions appear from the remaining terms, with coupling constants 
given by 



a?* = coaOwaf^ 



t ZR 

13 



1 BR 

13 



Pl R = C0S e w p B 

3-B „;„ oW 



ftf R = -em9wp% 



-y « - VY l^lj 

$l d L = cos 6 W 0? - sin 8 W ftW /3g L = - sin 8 W (3* - cos 6 W p% . (14) 



It is a simple matter to obtain the Feynman rules for the 7 lh h interactions from the lagrangian (they are 
identical in form to those obtained for the flavour- violating interactions gtc and g t u in refs. [Ij,ll5|). 
In figure [1] we present the Feynman diagrams for the decay [i — > e^f (in fact, for any decay of the type 






FIG. 1: (i — > e*y with effective anomalous vertices involving the couplings a, ft and 5. 



hi — > hi) with vertices containing the effective couplings a, ft and 8. Interestingly, the 8 contributions 
cancel out, already at the level of the amplitude [16J. The calculation of the remaining diagram is quite 
simple and gives us the following expression for the width of the anomalous decay lh — > hi in terms of 
the a and ft couplings: 



r(/ h - hi) 



,,7-R 



m hU a ih + a hi 



■yR* 1 2 



\ a lh + a hf*\ 2 ) 



64ttA 4 

+ 8m hv Im(<f - aj^ftl* - aJ h L ft? h * + afifi[ h ) 



7 |2 



W) 



(15) 



So, for the decay \x — > e7, using the data from [5j, we get (with A expressed in TeV) 

0.22 



BR(/i -> &i) = 



A 4 



(K + <T + K + <*l 2 ) + 4.3 x io 7 (i^i 2 + \^f) 



+ 1.3 x 10 4 Im« ftl e - a£P£ - a^ftj* + < /?!) 



(16) 



Now, all decays lh — * hi are severely constrained by experiment [B| especially in the case of /1 — > ei but 
also in t — ► ci and r — > fj/y. To obtain a crude constraint on the couplings, we can use the experimental 
constraint BR(^ — > e7) < 1.2 x 10" 11 and set all couplings but one to zero. With this procedure we get 
the approximate bound 



*efi I 



< 7.4 x 10" 



TeV" 



(17) 



and identical bounds for the oC 1 ^^ couplings. The constraints on the ft constants are roughly four orders 
of magnitude smaller. Using the same procedure for the two remaining LFV processes we get 



jL,R\ 

er i 



1.6 x 10~ 3 TeV" 2 



x lL,R\ 

flT I 

"A2 - 



< 1.3 x 10" 



TeV" 



(18) 



G 



with the P couplings even more constrained in their values. 

The experimental bounds on the various branching ratios are so stringent that they pretty much curtail 
any possibility of these anomalous operators having observable effects on any experiences performed at 
the 1LC. To see this, let us consider the flavour- violating reaction 77 — > Ihh, which in principle could 
occur at the ILC @. There are five Feynman diagrams involving the {a , (3} couplings that contribute to 
this process. There are also three diagrams involving the 5 couplings of eq. ©, but their contributions 
(once again) cancel at the level of the amplitude. The calculation of the cross section for this process 
is laborious but unremarkable. The end result, however, is extremely interesting. The cross section is 
found to be 

= r(i fc -i l7 ) , (19) 

dt rrih 3 s (rrih 2 — t) t (rn^ 2 — u) u 

with a function F 77 given by 

F 77 = m h ° (t + u)-12mltu + ml (t + u) (i 2 + 13 1 u + u 2 ) - m\ t u (7 t 2 + 24 1 u + 7 u 2 ) 

+ 12m 2 h t 2 u 2 (t + u)-6t 3 u 3 . (20) 

The remarkable thing about eq. l[T9|) is the proportionality of the (differential) cross section to the width of 
the anomalous decay 1^ — > li 7, which is to say (modulus the total width of //,, which is well known), to its 
branching ratio. A similar result had been obtained for gluonic flavour-changing vertices in refs. Il5j. 
Because the allowed branching ratios for the lh — ► li 7 are so constrained, the predicted cross sections 
for the ILC are extremely small. We have 

(7(77 -> fjT e+) ~ 1CT 8 x BR(^i -> ej) pb 
ct(77 — > t~ fi + ) — 10~ 5 x BR(r — * /17) pb 

(7(77 -> T~e + ) - 1CT 5 x BR(r ->■ e 7 ) pb , (21) 

with y/s = 1 TeV. With the current branching ratios of the order of 10 -12 for the muon decay and 10 -7 
for the tau ones, it becomes obvious that these reactions would have unobservable cross sections. 



Our conclusion is thus that the a?- and couplings are too small to produce observable signals in 
foreseeable collider experiments. However, both {otjj , 0Jj} and {afj , (3fj} are written in terms of the 
original {afj' W , ftfj' W } couplings, via coefficients (sine and cosine of 6\y) of order 1. Hence, unless there 
was some bizarre unnatural cancellation, the couplings {7 , Z} and {B , W} should be of the same order 
of magnitude. Since we have no reason to assume such a cancellation, we come to the conclusion that 
the a and (3 couplings are simply too small to be considered interesting. They will have no bearing 
whatsoever on anomalous LFV interactions mediated by the Z boson. From now on, we will simply 
consider them to be zero, which means that there will not be any anomalous vertices of the form 7 U lj. 



B. A set of free parameters 



In the previous section we have presented the complete set of operators that give contributions to the 
flavour violating processes e + e~ — » Ihh- However, these operators are not, a priori, all independent. 
It can be shown that (see refs. [n], [13, [I3> Ell f° r details), for instance, there is a relation between 
operators of the types O e B4> and O e B and some of the four fermion operators, modulo a total derivative. 
These relations between operators appear when one uses the fermionic equations of motion, along with 
integration by parts. They could be used to discard operators whose coupling constants are a and /3, 
or some of the four fermion operators. We used this argument to present the results in refs. [Til Il5l . [l8| 
in a more simplified fashion. However, in the present circumstances, we already discarded the a and 
(3 operators due to the size of their contributions to physical processes being extremely limited by the 
existing bounds on flavour-violating leptonic decays with a photon. Since we already threw away these 
two sets of operators, we are not entitled to use the equations of motion to attempt to eliminate another. 
Notice also that in most of the work that was done with the effective lagrangian approach one replaces, 
at the level of the amplitude, operators of the type Od c by operators of the type O e z^> by using Gordon 
identities. In fact, it can be shown that the following relation holds for free fermionic fields, 

el 9" ei = mj el ^ 4 - A a» a d a e? R . (22) 
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Notice that the use of Gordon identities is not the same thing as using the field's equations of motion to 

eliminate operators: in the latter case, one proves that different operators are related to one another and 

use those conditions to choose among them; in the former, all we are doing is re-writing the amplitude 

in a different form. And in our case, this procedure does not bring any simplification. 

Finally, using the equations of motion, a relation can be established between operators Od c and 0£> e , 

namely 

De +0 Be + (£ L e R ) \rle R t L +T u q L eu R + T\d R q L ^ =0 , (23) 

where the T e coefficients are the leptonic Yukawa couplings and e the bidimensional Levi-Civita tensor. 
We see that the relationship between these two operators involves four-fermion terms as well. This 
relation means we can choose between one of the two operators Ou s and O q , given that the four- 
fermion operators appearing in this expression have already been considered by us. This means that only 
one of the rjf- and vfc- couplings will appear in the calculation. We chose the first one and will drop, from 
this point onwards, the superscript "i?". Also, after expanding the operators of eq. ((3]), we see that the 
couplings always appear in the same combinations. We therefore define two new couplings, Or and 6l, 
as 

6 R = 0* + C , L = 0f « + Off* 9^ 9ff> . (24) 

As an aside, we must add that the use of equations of motion to simplify the effective lagrangian is not 
followed by all authors. For instance, the authors of [19( do not use them and consider instead a fully 
general set of dimension six effective operators. 
The original lagrangian is now reduced to 

C = C eeh i h + Dc + O 0e + 0$ + 0$ + + h.c. . (25) 

We will not present the Feynman rules for the four fermion interactions because they are obvious and 
rather cumbersome to write. The remaining Feynman rules we will use are presented in figure [2] where 
lepton momenta follow the arrows and vector boson momentum is incoming. For completeness, we 
included the afj and in this figure, but we remind the reader that we have set them to zero. 



J2 [%7r(um R P» + a fh R '^) + 7„7L(af,V + uf h L " q v ) 
K a WV^ + v <v {(ilan + PH'jl)] (^ g va - k v g» a )+ 

J2 [(lR Vhl Pa + 7L nth la) + V J a (6 L J L + 6 R Jr)} 




A-' 



FIG. 2: Feynman rules for anomalous Z lh U and U lh vertices. 



IV. DECAY WIDTHS 

As we said before, all LFV processes are severely constrained by experimental data. Now that we have 
settled on a set of anomalous effective operators, we should first consider what is the effect of those 
operators on LFV decays. The existing data severely constrains two types of decay: a heavy lepton 
decaying into three light ones, lh — > III, such as t~ — * e~e + e~, and decays of the Z boson to two 
different leptons, Z — > Ihh (such as Z — > t + e~). Flavour- violating processes involving neutrinos 
in the final state (such as, say, Z — > v r D e ) are not constrained by experimental data, as they are 
indistinguishable from the "normal" processes. 
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For the 3-lepton decay, there are three distinct contributions, whose Feynman diagrams are shown in 
figure [3] for the particular case of fi~~ — ► e _ e + e~. As before, the contributions involving the S operators 
cancel at the level of the amplitude and have absolutely no effect on the physics. Using the Feynman 
rules in figure [2] and the four fermion lagrangian we can determine the expression for the decay lh —* III. 
Remember that I stands for a massless lepton whatever its flavour is. The decay width obtained is the 




FIG. 3: Feynman diagrams for the decay \x 



sum of three terms, to wit 

T(l h -> III) 



r 4f {i h -> in) + r z (i h -> in) + r mt (i h — > in) 



(26) 



where T^f contains the contributions from the four- fermion graph in figure El T z those from the Feynman 
diagram with a Z boson and Ti nt the interference between both diagrams. A simple calculation yields 



where 



Tz(k 



ill) 
III) 



6144 tt 3 A 4 
768M 4 tt 3 A 4 



[\S LR \ 2 + \S RL \ 2 +4(\V LL \ 2 + \V RR \ 2 )] 



+ 



in 



10M? 



(I^lI 2 + IN 2 ) v 2 m\ + -Re ( mh 9* L + i lhl R )vm 6 h 



|^| 2 )M Z 2 + 6(|^| 2 + |^| 2 )« 2 1 



Fint(lh —> HI) 



v ml 



768M z 2 tt 3 A 4 

(gv 



3m| 
lOAf? 



(g v + g A )Re(e L Vl L ) + 



m h 



g A )B e (9 R V* R ))--J±[l 



{gv - 9 a) Re {r] M V LL ) 



(gv + 9a) Re (mhV RR ) + 



gv 



sin 9w cosdyy 



\ + sin^ 



9A = 



K 

ism 9w cos9w 



(27) 



(28) 



and e is the elementary electric charge. An important remark about these results: they are not, in fact, 
the exact expressions for the decay widths. The full expressions for T z{hi —* III) and Ti nt (lh III) are 
actually the sum of a logarithmic term and a polynomial one. However, it so happens that the first four 
terms of the Taylor expansion in m,h/M z of the logarithm cancel the polynomial exactly. The expressions 
of eq. ([2T| are therefore the first surviving terms of that Taylor expansion, and constitute an excellent 
approximation to the exact result, and one that is much easier to deal with numerically (the cancellation 
mentioned poses a real problem in numerical calculations). 
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As for the LFV decays of the Z-boson, there is an extensive literature on this subject @. There are, 
of course, no four fcrmion contributions to this decay width, and a simple calculation provides us the 
following expression: 

r(Z ^ Mi) = [(Mf-mD 2 (|w| 2 + |^| 2 )+4(^ + 2M z V(|^| 2 + |^| 2 ) 

+ 4m h (ml-M*)vRe(6 L r) h i + 6 R r)t h )] . (29) 



V. CROSS SECTIONS 



In this section we will present expressions for the cross sections of various LFV processes that may 
occur at the ILC. There are three such processes, namely (1) e + e~ — > fi~e + , (2) e + e~ — > r~e + and 
(3) e + e~ — > t~/j + , as well as the respective charge-conjugates. We have calculated all cross sections 
keeping both final state masses. However, given the energies involved, the contributions to the cross 
sections which arise from the lepton masses are extremely small, and setting them to zero is an excellent 
approximation. We thus present all formulae with zero leptonic masses, as they are much simpler than 
the complete expressions. In figures U and [5] we present all diagrams that contribute to the process 
e + e~ — > /i~e + . A brief word about our conventions. There are two types of LFV production cross 




FIG. 4: Feynman diagrams describing the process e e — > fi e 

sections, corresponding to different sets of Feynman diagrams. In the case of process (1), we see from 
figures [4] and [5] that the reaction can proceed through both a i-channel and an s-channel - this is obvious 




FIG. 5: Interpretation of the four fermion terms contributing to the process e e — » fi e in terms of currents; 
notice the analog of a t channel and an s one. 



for the diagrams involving the exchange of a photon or a Z boson. For the four-fermion channels less 
so, but figure [5] illustrates the t and s-channel analogy. Depending on the "location" of the incoming 
electron spinor in the operators of eq. I|5|). we can interpret those operators as two fermionic currents 
interacting with one another, that interaction is obviously analog to the two different channels. Process 
(2) has diagrams identical to those of process (1). Process (3), however, can only occur through the s 
channel - that is obvious once one realizes that for process (3) there is no positron in the final state. In 
fact for process (3) there are only "s-channel" contributions from the four-fermion operators. 

A simple way of condensing the different four-fermion cross sections into a single expression is to adopt 
the following convention: we will include indices "s" and "t" in the four fermion couplings. If we are 
interested in the cross sections for processes (1) and (2) - which occur through both s and t channels - 
then all "s" couplings will be equal to the "t" ones. If we wish to obtain the cross section for process (3) 
(which only has s channels) we must simply set all couplings with a "t" index to zero. We have further 
considered the likely possibility that in the ILC one may be able to polarize the beams of incoming 
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electrons and positrons [20|. Thus, 07,7 represents the polarized cross section for an I polarized electron 
and a J polarized positron, with {7 , J} = {R , L} - that is, beams with a right-handed polarization or 
a left-handed one. The explicit expressions for the four-fermion differential cross sections are then given 

by 



da LL 


1 


dt 


16tts 2 A 4 


d(JRR 


1 


dt 


16tts 2 A 4 


da LR 


1 


dt 


16tts 2 A 4 


ddRL 


1 


dt 


16tts 2 A 4 



[^ 2 \V£ L - 



V, 



v, 



LL I " 

t |2 
RR\ 



i 2 (4|^ 



LR\ 



Hi I 



LL 



\u 2 \S 



RR 



sl L \ 2 



sU 2 



-s 2 mi R \ 2 



\S RL \ 2 )] 



\sIr\ 2 ) 



\SIr\ 2 )] 



\S RL \ 2 )] 



(30) 



See appendix [B] for the full calculation. The unpolarized cross section is obviously the averaged sum over 
the four terms of eq. (|30"1) . To re-emphasize, the four-fermion cross section for processes (1) and (2) is 
obtained from this expression by setting all "s" couplings equal to the "t" ones; and to obtain the cross 
section for process (3) one must simply set all "t" couplings to zero. 
The total cross sections for each of the processes are then given by 



r (l,2) 



(e-e+ -> l h e+) 



r ^ 



(1,2) 
Z 

(3) , 



4/ 
T (3) , 



r (l,2) 



.(3) 



(31) 



where az is the cross section involving only the anomalous Z interactions of figure [2 174/ the four-fermion 
cross section - whose calculation we already explained - and <Tj nt the interference between both of these. 
The S couplings also present in figure [5] end up not contributing at all to the physical cross sections, once 
again. For completeness, then, the remaining terms in the differential cross section for processes (1) and 
(2) are given by 



da 



(1,2) 



dt 

with 

F\{gA,9v) 



v 2 [^( gA , gy )|g L | 2 + F 1 (g A ,-g v )\0 R \ 2 ] + F 2 (g A ,g v )\ Vlh \ 2 + F 2 (g A ,-g v )\r ]hl \ 

32 7r A 4 (M 2 - sf s 2 {M 2 - tf 



= 2 {{gA + gvf [st(2M* + 2uM 2 + s 2 + t 2 )~uM 2 (u M 2 + 2 s 2 + 2 1 2 )} + 
2 {g\ + g 2 v )u {t 3 + s 3 ) + (g A - g v ) 2 su 2 t) 



(32) 



s 2 + t 2 + s t) + 2 g A g v (M 2 - s) (M 2 - t)] . (33) 



F 2 (gA,gv) =-tus [(g 2 A +g 2 v )(3M* + 3uM : 
The interference term is given by 

dot? {t-s)v 2 

dt 16 tf A 4 (M 2 - s) s 2 {M 2 - t) 

(g A Re (9 L S* LR ~ 9 R S* RL ) + g v Re (9 L S* LR + 6 R S* RL )) (s t + (-M 2 + s + t)u) 
+ 4(g A - g v ) (s + t)u Re (9 L V£ L ) -4(g A + gv) (s+t)u Re (9 R V^ R ) 
For process (3), we have 



(34) 



daf -v 2 
dt ~ 32ttA 4 (M 2 2 - sf s 2 



2v 



[(g 2 A +g 2 v)(2tu-s 2 )(\9 L \ 2 + \9 R \ 2 ) 



+ 2g A9v s(t - u) (\9 L \ 2 - \9 R \ 2 )] - (g 2 A + g 2 v ) t u s (\ Vhl \ 2 + \ Vlh \ 2 ) 
and finally, the interference terms are 



dt 8ttA 4 (M 2 -s) s 2 



(gv - g A ) Re (9 L V£ L ) + (g v + g A ) Re (9 R V£ R ) 



(35) 



(36) 
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At this point we must remark on the different energy behavior that these various terms follow. Once 
integrated in t, the four-fermion terms grow linearly with s, whereas those arising from the anomalous 
Z couplings have a much smoother evolution with s - whereas the first ones diverge as s — > oo, the 
second ones tend to zero. See appendix IBI for the expressions of the integrated cross sections. This could 
be interpreted as a clear dominance of the four-fermion terms over the remaining anomalous couplings. 
However, we must remember that we are working in a non-renormalizable formalism. We know, from the 
beginning, that these operators only offer a reasonable description of high-energy physics up to a given 
scale, of the order of A. The dominance of the four-fermion cross section must therefore be carefully 
considered - it may simply happen, as there is nothing preventing it, that the four-fermion couplings of 
eq. ^ are much smaller in size than the Z boson ones of figure^ 

As we saw, the 6 couplings end up not contributing to either decay widths or cross sections (and 
this is true regardless of whether the light leptons are considered massless or not). As we mentioned 
before, their inclusion could be interpreted as an on-shell renormalization of the leptonic propagators. 
On that light, their cancellation suggests that the effective operator formalism is equivalent to an on- 
shell renormalization scheme. This is further supported by the fact that the list of effective operators 
of ref. [llj was obtained by using the fields' equations of motion to simplify several terms. However, 
we must mention that at least in some Feynman diagrams (some of those contributing to 7 7 — > luhi 
for instance), the "<5-insertions" were made in internal fermionic lines, so that this cancellation is not 
altogether obvious. 



A. Asymmetries 



In a collider with polarized beams, asymmetries can play a major role in the determination of flavour- 
violating couplings. A great advantage of using these observables is that, as will soon become obvious, 
all dependence on the scale of unknown physics, A, vanishes due to their definition. There is a strong 
possibility that the ILC could have both beams polarized, therefore the measurement of polarization 
asymmetries could be very interesting. For a more detailed study see (20| . A particulary appealing 
situation is found when the contributions from the Z boson anomalous couplings are not significant 
when compared with the four fermion ones. In this case the study of asymmetries would allow us, in 
principle, to determine each four-fermion coupling individually. We will now concentrate on one of the 
most feasible scenarios, which is to have a polarized electron beam and an unpolarized positron beam. We 
will take both the right-handed and left-handed polarizations to be 100%, which is obviously above what 
is expected to occur (recent studies show that a 90 % polarization is attainable) [2(J. The differential 
cross sections for left-handed (P e - = —1) and right-handed (P e - = +1) polarized electrons are 



{4u*\Vi L + V* L \ 2 +1?\S t RL \ 2 + * l \S> L 

(4 u 2 |^ + ^«| 2 + ^ 2 |5i«| 2 + S 2 |^ L | 2 ) • (37) 



dt 16tts 2 A 4 



Two forward-backward asymmetries for the left-handed and right-handed polarized cross sections can 
now be defined as 

Jo /2 da L(R) (9) - /; /2 da L(R) (9) 

A FB,L(R) = V 6 °) 

and we can also define a left-right asymmetry, given by 

A LR = ^-^ , (39) 
cl + ctr 

where <Jl(r) is the total cross section for a left-handed (right-handed) polarized electron beam. Note that 
we have assumed that the polarization of the final state particles is not measured. Otherwise we could 
get even more information by building an asymmetry related to the measured final state polarizations. 
Using the expressions on appendix [51 it is simple to find, for these asymmetries, the following expressions: 

A = V\V£l + VIl\ 2 -3\S rl \ 2 

FB ' L 16|V/ L + V2J 2 +4|^ L P + 12|S£ fl P 1 > 
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and 



Apr r = '"^ 1 li^i (41) 



Finally, the left-right asymmetry reads 

, = \V£l + vL\ 2 ~ \Y£r + vk R \ 2 + l^kl 2 - l^l 2 + 3 (|S£ fl | 2 - \S RL \* 
LR \V£ L + V[ L V + \Vj> R + V RR \* + l^kl 2 + l^l 2 + 3 (\S iR \* + \S RL \* 



which has no dependence on A. Notice that all of these expressions assume an unpolarized positron 
beam, and a completely polarized electron beam, either left- or right-handed. If the electron beam is not 
perfectly polarized, but instead has a percentage of polarization P e - , we can still write 

ap e _ = <t [1 ~P e -A LR ] (43) 

with Co = {?l + 0\r)/4. So if in reality we only have access at the ILC to beams with +80 % (- 80 %) 
polarization we could still use them to determine ctq and Al R . If we had access to a positron polarized 
beam, we could then write a similar expression for the cross section obtained from the polarized positrons. 
Notice that Alr would be different - the indices left and right would then refer to the positron and not 
to the electron. 

The most interesting possibility is, of course, when both beams are polarized, with different percent- 
ages, P e - and P e + . We could then perform experiments where the four different combinations of beam 
polarizations were used. The resulting cross section would be 
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(l + P e -)(l + P e+ )a RR + {l-P e -)(l-P e+ )a LL + 



(1 + P e -)(1 - P e+ )a RL + (l-P e -)(l + P e+ )a LR 



(44) 



As such, we would be able to determine o RR , <tll, &rl and ulr - and consequently each of the four 
four-fermion couplings, V LLl V RR , S RL and S LR . 



VI. RESULTS AND DISCUSSION 



In the previous sections we computed cross sections and decay widths for several flavour-violating pro- 
cesses. We will now consider the possibility of their observation at the ILC. To do so we will use one set 
of parameters [2(| for the ILC, i.e., a center-of mass energy of i/s = 1 TeV and an integrated luminosity 
of C = 1 ab _1 . At this point we remark that, other than the experimental constraints on the flavour- 
violating decay widths computed in sec. UVK see tabled]), we have no bounds on the values of the anomalous 
couplings. The range of values chosen for each of the coupling constants was 10 -4 < |a/A 2 | < 10 _1 , where 
a stands for a generic coupling and A is in TeV. For a w 1 the scale of new physics can be as large as 
100 TeV. This means that if the scale for LFV is much larger than 100 TeV, it will not be probed at the 
ILC unless the values of coupling constants are unusually large. The asymmetry plots are not affected 
by this choice as explained before. We will therefore generate random values for all anomalous couplings 
(four-fermion and Z alike) , and discard those combinations of values of the couplings for which the sev- 
eral branching ratios we computed earlier are larger than the corresponding experimental upper bounds 
from table fl] This procedure allows for the possibility that one set of anomalous couplings (the Z or 
four-fermion ones) might be much larger than the other. When an acceptable combination of values is 
found, it is used in expressions (|3H| - (f3"5| to compute the value of the flavour-violating cross section. In 
figure [S] we plot the number of events expected at the ILC for the process e + e~ — > r~ e + , in terms of 
the branching ratio BR(Z — > re). To obtain the points shown in this graph, we demanded that the 
values of the effective couplings were such that all of the branching ratios for the decays of the r lepton 
into three light leptons and BR(Z — ► r fx) were smaller than the experimental upper bounds on those 
quantities shown in table |TJ We observe that, even for fairly small values of the r flavour-violating decay 
branching ratios (10 -9 -10 -6 ), there is the possibility of a large number of events for the anomalous cross 
section. 

By following the opposite procedure - requiring first that the branching ratios BR(Z — > re) and 
BR(Z — > t fi) be according to the experimental values, and letting BR(r — > III) free, where I is 
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TABLE I: Experimental constraints on flavour- violating decay branching ratios 

Process Upper bound 

2.0 x 1(T 7 



r — > eee 
t — > e [i [i 

T 



■7 



2.0 x 10" 

^tee 1.1 x 10~ 7 

t — > fj,fifi 1.9 x 10 -7 

H -> eee 1.0 x 10 -12 





— > e /i 


1.7 x 10 -6 




^ CT 


9.8 x 10~ 6 




-> T/i 


1.2 x 10" 5 




either an electron or a muon - we obtain the plot shown in figure [Jj This time we analyse the process 
e + e~ — * t~ but a similar plot is found for e + e _ — > t~ e + . The number of events rises sharply with 
increasing branching ratio of r into three leptons. It is possible to discern a thin "band" of events in the 
middle of the points of figure [71 rising linearly with BR(r — * III). This "band" corresponds to events 
for which the four-fermion couplings are dominant over Z° events. In that case, they dominate both 
BR(t — > III) and a(e + e~ — ► r~ fi + ), and the larger one is, the larger the other will be - which explains 
the linear growth of this subset of points in the plot of figure [JJ This "isolated" contribution from the 
four-fermion terms is not visible in figure [S] since the branching ratios of the Z decays are independent 
of those same couplings. Finally and for completeness, in figure [8] we show the values of the asymmetry 
coefficient Afb.r defined in (j¥T|) . for the process e + e~ — > t~ versus the three-lepton decay of the 
r. A similar plot is obtained for the asymmetry Afb,l- We observe a fairly uniform dependence on the 
branching ratio BR(t — > III), which is to say, on the values of the four-fermion couplings. However, 
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BR(t -> 1 1 1) 



FIG. 7: Number of expected events at the ILC for the reaction e e — > r fj, , with a center-of-mass energy of 
1 TeV and a total luminosity of 1 ab" 1 . 
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Afb,r asymmetry for the process e + e — > r /i + versus BR(t III) 
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there is a significant concentration of "points" near the maximum allowed value for this cross section, 
0.75. 

Finally, we also considered another possible process of LFV, namely 7e~ — > fi~ Z°. There are three 
Feynman diagrams contributing to this process, one of which involving a quartic vertex which emerges 
from the effective operators of eqs. © and This process might occur at the ILC, if we consider 
the almost-collinear photons emitted by the colliding leptons, well described by the so-called equivalent 
photon approximation (EPA) [2l[. An estimate of the cross section for this process, however, showed it to 
be much lower than the remaining ones we considered in this paper. This is due to the EPA introducing 
an extra electromagnetic coupling constant into the cross section, and also to the fact that the final state 
of this process includes at least three particles (one of the beam particles "survives" the interaction)- 
thus there is, compared to the other processes which have only leptons in the final state, an additional 
phase space suppression. Notice, however, that an optional upgrade for the ILC is to have e 7 collisions, 
with center-of-mass energies and luminosities similar to those of the e + e~ mode, so this cross section 
might become important. 

The flavour-violating channels are experimentally interesting, as they present a final state with an 
extremely clear signal, which can be easily identified. The argument is that the final state will always 
present two very energetic leptons of different flavour, more to the point, an electron and a muon. LFV 
can be seen in one of the three channels e + e~ — > fjT e + , e + e~ — > r~e + , e + e~ — > r~/i + and charge 
conjugate channels. The first channel is the best one, with the two leptons back to back and almost free 
of backgrounds. For the other production processes, we may "select" the decays of the tau that best 
suit our purposes: for the second we should take the tau decay t~ — > fi~i7^i> Tl and for the third process, 
t~ — ► e~v e v r . The branching ratios for both of these tau decays are around 17%, so the loss of signal is 
affordable. The conclusion is that, for every lepton flavour-violating process, one can always end up with 
a final state with an electron and a muon. If the ILC detectors have superb detection performances for 
these particles, then the odds of observing violation of the leptonic number at the ILC, if those processes 
do exist, seem reasonable. 

Clearly, our prediction that significant numbers of anomalous events may be produced at the ILC needs 
to be further investigated including the effects of a real detector. Notice also that due to beamstrahlung 
effects which reduce the effective beam energy, the total LFV event rates might be reduced, specially 
in the case of the four-fermion cross sections, which increase with s. Also, one must take into account 
the many different backgrounds that could mask our signal. And the fact that, even in the best-case- 
scenario, only a few thousand events are produced with an integrated luminosity of 1 ab _1 , could limit the 
signal-to-background ratio. A careful study of the background to LFV processes lies clearly beyond the 
scope of this paper. Nevertheless, we will show that with some very simple cuts most of the background 
can be eliminated. Due to the weaker experimental constraints on processes involving r leptons, the 
most promising LFV reactions at the linear collider are fj,e and fir production. For illustrative purposes 
we will study the backgrounds to the LFV process e + e~ — > r + e~ — > /i + e~ v T . The main sources 
of background to this process are e + e~ — > /i + e~ v e and e + e~ — > t + t~ — > /i + e~ v e v T v T . The 
cross section to the background process e + e _ — » /i + e~ v e was calculated using WPHACT [22j and 
confirmed using Racoon WW [23j. The cross section for the remaining background was evaluated using 
PYTHIA [24]. In e + e~ — > t + e~ — > /i + e~ v T the electron is produced in a two body final state. 
Therefore its energy is approximately half of the center of mass energy. Furthermore, if 9 e is the angle 
between the electron and the beam, then the transverse momentum of the electron is px — sin(0 e ). 
This means that a cut in 9 e implies a cut in px- The main contribution to this cross section comes from 
the four-fermion interaction. There are no propagators involved and consequently the dependence in 9 e 
(and in px) is very mild. This can be seen from the expression (|B5p in the appendix. Making all coupling 
constants Vij and Sij equal, it can be shown that a 10 degree cut will reduce the cross section by 2 % 
while a 60 degree cut will reduce it only by 58 %. In table I we show the cross sections for the signal 
and for the backgrounds as a function of a cut in 9 e and a corresponding cut in px- For the signal we 
start with a cross section of 5 fbarn when no cuts are applied. Due to the mild dependence on 9 e , a cut 
of 60 degrees will make the signal well above background. A further cut on the energy of the electron 
could be applied, say E e > 300 GeV. This would not affect the signal but will reduce the background 
even further. All calculations were performed at tree level with Initial State Radiation and Final State 
Radiation turned off. Another possibility for background reduction would be to use the polarisation of 
the beams, a method known to be very efficient. Notice, however, that this procedure might affect the 
extraction of four-fermion couplings from polarised beam experiments - if the signal is observed only for 
certain combinations of beam polarisations, it could happen that only certain couplings, or combinations 
thereof, can be measured. 
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Cut in 6 e (degrees) 10 20 30 40 50 60 



e + e" -> t+ e~ ->■ fj,+ e~ u„ u T 4.9 4.6 4.1 3.5 2.8 2.1 

e+e" -> fi + e" u^u E 68.2 26.3 10.8 4.4 1.6 0.5 

e+e- -> r+ r" -» /z + e~ z/ M ^ e i/ T i/ T 1.3 0.8 0.3 0.2 0.06 0.01 



TABLE II: Cross sections (in fbarn) for the LFV signal and most relevant backgrounds to that process for several 
values of the angle cut between the the outgoing electron and the beam axis. 



Finally, some comments on the dependence of these results vis-a-vis expected improvements on the 
measurements of the LFV branching ratios of table [J Could it be that future experiments would tighten 
the constraints so much that there was no room available for discovery? Tau physics at BABAR and 
BELLE has provided the best limits so far on LFV involving the r lepton. The combined results from 
BABAR and BELLE on r — > Ij are now reaching the level of 10 -8 and will be close to just a few 10~ 8 
by 2008 [2|| . More important to us are the decays r — > III, due to the constraints imposed on the four 
fermion operators. The latest results on Br[r — > III) from BABAR and BELLE are of the order of 10~ 7 , 
with less than 100 /6 _1 of data analysed. A value of the order of a few 10 -8 is expected when all data is 
taken into account [26| . Other planned experiments like MEG or SINDRUM2 (see [27]] ) will provide much 
more precise results for both /i — > and /i e conversion, respectively. However, those results will not 
constrain any further the four-fermion couplings. The current limit Br(fi -teee) < 10~ 12 at 90% CL f29j 
already excludes the possibility of finding LFV in the \ieee coupling. This limit will be improved by 
the Sundrum experiment (see [27|]). Another possibility is the GigaZ option for the ILC, which probably 
would be earlier than an energy upgrade to 1 TeV. Again, the limits on the LFV branching ratios of 
the Z boson would be improved [28| but the bounds on the four-fermion couplings would not be affected. 
Lastly, LFV searches will also take place at the LHC. Preparatory studies on the LFV decay r — > ^ M M 
are being conducted by CMS [3(| , ATLAS [3lJ and also by LHCb [32| . During the initial low luminosity 
runs (10-30 fb~ 1 /ye&i) for 2008-2009, searches for this decay may be possible. So far the limits predicted 
are only slightly better than the known limits from the B-factories. Therefore, in the foreseeable future, 
the constraints on the four fermion r couplings arising from the branching ratios of table [J could go down 
one order of magnitude, to be of the order of 10 -8 . Accordingly, and repeating the calculations that led 
to figs. [6] and the maximum number of events expected at the ILC also goes down by one order of 
magnitude, to about 1000 events. Given the discussion on backgrounds above, we expect that detection 
of LFV at the ILC would still be possible, although harder. 
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APPENDIX A: SINGLE TOP PRODUCTION VIA GAMMA-GAMMA COLLISIONS 



In section fill Al we argued that the couplings corresponding to the operators of eqs. (U) were extremely 
limited in size by the existing experimental data for the branching ratios of the decays lh — > h 7- In fact, 
we even showed that the cross sections for the processes 77 — > lh k, eq. (p~9|) , were directly proportional 
to those branching ratios, and their values at the ILC were predicted to be exceedingly small. It is easy 
to understand, though, that we can define operators analogous to those of eqs. ^ for quarks instead of 
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lcptons. In particular, we can consider flavour-changing operators involving the top quark, which would 
describe decays such as t — ► U7 or ( — * 07 - and these decay widths have not yet been measured. 
More importa ntly , their values may vary immensely, depending on the model one uses to calculate them. 
According to [33], the branching ratios for these decays range from their SM value of ~ 10~ 16 (for the u 
quark), ~ 10™ 12 (for the c quark) to ~ 10~ 6 (for both quarks) in supersymmetric models with R-parity 
violation. The total top quark width being also a lot larger than the tau's or the muon's, it seems possible 
that the cross section for single top production via flavour-violating photon-photon interactions presents 
us with observable values. 

The corresponding calculation is altogether identical to the one we presented for the leptonic case. We 
find an expression for the width of the anomalous decay t — > q 7 similar to that of eq. (|15[) , 



r(t -> qi ) 



m?(|< + <*| 2 + \*f + <T) + 16« 2 m t \ 2 + K\ 2 ) 

(Al) 



64ttA 4 

+ 8 m t v Im(a]*01 - al^Pj* - alt + ®tqPta) 



with new couplings a and f3 (we re-emphasize that these new couplings are not in the least constrained 
by the arguments we used in section fill Ap and q = {u , c}. Likewise, considering that the top quark's 
charge is 2/3 and the quarks have three colour degrees of freedom, we may rewrite the analog of eq. (fll?)) 

as 

uo Jraf s (m t — t) t (nit — u) u 

with F-yy given by an expression identical to eq. I|20[). with the substitution vtih <-> mt. With a top total 
width of about 1.42 GeV and for y/s equal to 1 TcV, this expression can be integrated in t (with a pt cut 
of 10 GeV on the final state particles, to prevent any collinear singularities) and the total cross section 
estimated to be of the order 

ct(77 -> tq) - 90 x BR(i -> 97) pb . (A3) 

We see a considerable difference vis-a-vis the predicted leptonic cross sections, from eqs. (|2"Tj) - this one is 
much larger. To pass from the photon-photon cross section to an electron-positron process, we apply the 



standard procedure: use the equivalent photon approximation 21( to provide us with the probability of 
an electron/positron with energy E radiating photons with a fraction x of E and integrate eq. (|A2[) over 
x. For recent studies of photon-photon collisions at the ILC, see for instance (34|. The numerical result 
we found for the single top production cross section is 

cr(e + e~ -> e+ tq) = 1.08 x BR(i -> qj) pb . (A4) 

For an integrated luminosity of about 1 ab _1 , this gives us about one event observed at the ILC for 
branching ratios of t — > qj near the maximum of its theoretical predictions 35], ~ 10~ 6 . Clearly, this 
result means that this process should not be observed at the ILC, even in the best case scenario. However, 
in the event of non-observation, eq. (|A4[) could be useful to impose an indirect limit on the branching 
ratio BR(£ — > (77). Several authors have studied single top production in e + collisions [I2T [36|. For 
gamma-gamma reactions, single t op p roduction at the ILC in the framework of the effective operator 
formalism may has been studied in j37| . and for specific models, such as SUSY and technicolor, in ref. (34|. 



APPENDIX B: TOTAL CROSS SECTION EXPRESSIONS 



We write the amplitude for the four fermion cross sections in two parts. One for the s channel and the 
other one for the t channel. In doing so we are generalizing the four fermion lagrangian which for a gauge 
theory has equal couplings for both s and t channels. For the s channel the amplitude reads 

T ij = ji [V£ (v e iai i u e )(ui h 'y a 'yjv ll ) + (w e 7*«e)(uj h 7ji;j ( )] (Bl) 
while for the t channel we have 

= - -rg [Vfj (ui h j a %u e )(v e 'y a j : jv e ) + Sjj (ui h ^iU e ){v e ^jV e )] (B2) 
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with i,j = L, R. With these definitions we can write 
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(B3) 

and to obtain the expressions when only the t or s channels are present, you just have to set the s 
couplings or the t couplings, respectively, equal to zero, u, t and s are the Mandelstam variables defined 
in the usual way. 

The cross sections for polarized electron and positron beams with no detection of the polarization of 
the final state particles were given in eq. ([50]) . The International Linear Collider will have a definite 
degree of polarization that will depend on the final design of the machine. For longitudinally polarized 
beams the cross section can be written as 

d(Tp <- P < + ^(1 + Pe-)(1 + P e+ )^ + (l-P e -)(l-P e+ ) d<JLL ' 



dt 



dt 



(1+P e -)(1-P e+ ) 



d<JRL 
dt 



(1-P e -)(1 + P e+ ) 



dcr LR 
dt 



(B4) 



where <jrl corresponds to a cross section where the electron beam is completely right-handed polarized 
(P e - = +1) and the positron beam is completely left-handed polarized (P e + = —1). This reduces to 
the usual averaging over spins in the case of totally unpolarized beams. For the general expression for 
polarized beams, as well as a study on all the advantages of using those beams, see (20j . 

In the main text we presented expressions for the differential cross sections. For completeness we now 
present the formulae for the total cross sections. For the four-fermion case, the expressions have a very 
simple dependence on the px cut one might wish to apply, so we exhibit it. The quantity x = y/l — kp\ / s, 
with pt being the value of the minimum transverse momentum for the heaviest lepton, gives us an 
immediate way of obtaining these cross sections with a cut on the pr of the final particles. The total 
cross section is obviously the sum over all polarized ones, which gives us 
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(B5) 

As explained in the main text, the cross sections for processes (1, 2) are obtained from eq. (|B5[) by setting 
all of the "s" couplings equal to the "t" ones, and, for process (3), by setting the "t" couplings to zero. 

For the remaining cross section expressions we imposed no pt cut on any of the final particles. The 
total cross section for the Z couplings is given by, for processes (1, 2), 

„,2 
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with interference terms 

(1,2) = v 2 

a int - 48A 4 tt(s-M z 2 ). 

* (9A - 9v) {(12M 4 + 6sM 2 z - 14s 2 )Re (9 L V£ L ) - s 2 Rc (6 R S* RL )} 
- (9A + 9v) {(12M Z 4 + 6sM 2 z - Us 2 )Re (9 R V£ R ) - s 2 Rc (6 L S* LR )} 

{9 A - gv) { (4M 4 + 8sM 2 - 4s 2 )Re (6 L V£ L ) - s 2 Rc (6 R S* RL )} 
1 {(4M 4 + 8sM 2 - 4s 2 )Rc (9 R V^ R ) - s 2 Re (9 L S* LR )} 

Finally, for process (3), we have 
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APPENDIX C: NUMERICAL VALUES FOR DECAY WIDTHS AND CROSS SECTIONS 

We present here numerical values for the several decay widths and cross sections given in the text. We 
have set, in the following expressions, A equal to 1 TeV, the dependence in A being trivially recovered if 
we wish a different value for it. 
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For the cross sections, taking y/s = 1 TeV and imposing a cut of 10 GeV on the pr of the particles in 
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the final state, we have (in picobarn): 
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